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ABSTRACT 

We study the kinematics of gaseous disks in triaxial dark matter halos using the 
04 ' closed-loop orbit solutions in non-axisymmetric potentials. The orbits are in general 

non-circular and, for given triaxiality, their ellipticity depends on the ratio of escape 
to circular velocities, V^ sc /V^. This ratio increases steeply towards the center for cold 
dark matter (CDM) halo density profiles, implying that even minor deviations from 
spherical symmetry may induce large deviations from circular orbits in the velocity 
■ field of a gaseous disk, especially near the center. This result suggests that caution 

\q [ should be exercised when interpreting constraints on the presence of density cusps in 

. the dark halo derived from the innermost velocity profile. Simulated long-slit rotation 

curves vary greatly in shape, depending primarily on the viewing angle of the disk 
and on its orientation relative to the principal axes of the potential. "Solid-body" 
rotation curves — typically interpreted as a signature of a constant density core in the 
dark matter distribution — are often obtained when the slit samples velocities near 
the major axis of the closed loop orbits. Triaxial potentials imprint specific symme- 
tries in 2D velocity fields, generally inducing "twists" in the isovelocity contours and 
anti-symmetric patterns in opposite quadrants. We suggest that triaxial halos may be 
responsible for the variety of shapes of long-slit rotation curves of low surface bright- 
ness (LSB) galaxies, as well as for the complex central kinematics of LSBs, which are 
sometimes ascribed to the presence of "radial motions" in the gas. We argue that LSB 
rotation curves might be reconciled with the structure of CDM halos once the effects 
of halo triaxiality on the dynamics of gaseous disks are properly taken into account. 

Key words: cosmology: dark matter - galaxies: formation - galaxies: kinematics and 
dynamics 
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1 INTRODUCTION 

It is commonly believed that the inner regions of low surface 
brightness (LSB) galaxies are ideal probes of the inner struc- 
ture of dark matter halos. Although estimates of the stellar 
mass-to- light ratio remain somewhat uncertain, the baryonic 
component contribution to the mass budget in these galax- 
ies is generally thought to be small. Under this assumption, 
dynamical tracers of the potential such as rotation curves 
are expected to cleanly trace the dark matter distribution 
in LSBs. This provides important astrophysical clues to the 
nature of dark matter, since the spatial distribution of dark 
material in these highly non-linear regions is expected to be 
quite sensitive to the physical properties of the dark matter. 

LSB rotation curves may thus be contrasted directly 
with theoretical predictions of the inner structure of ha- 
los, and there is now an extensive body of work in 



the literature that reports substantial disagreement be- 
tween the shape of LSB rotation curves and the circu- 
lar velocity curve s of simulated cold dark matter (CDM) 
halos (see, e.g., |Flfjres_^njl^rirnaclJ 11994 iMoorel 11994 ; 
IMcGaugh and de BlokHl99d Ide Blok et alJl20oC Some of 
these rotation curves are fit better by circular velocity curves 
arising from density profiles with a constant density "core" 
rather than by the "cuspy" densi ty profiles commonly used 
to fit the structure of CDM halos iNavarro et alll99rj.ll997l . 
hereafter NFW). This discrepancy adds to a growing list of 
concerns regarding the consistency of CDM with observa- 
tional constraints on the scale of individual galaxies (see, 
e.g.. lSellwood and Kosowskvll200ll) that has prompted calls 
for a radi cal revision of the CDM para digm on small scales 
(see, e.g.. lSpergel and Steinhard3l2OO0() . 

Before accepting the need for radical modifications to 
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CDM it is important to note a number of caveats that apply 

to the LSB rotation curve problem. For instance, many of 

the early rotation curves where the disagreement was noted 

were significantly aff ected by beam smearing in the HI data 

iS waters et all boOO). The observational situation has now 

improved substantially thanks to higher-resolution rotation 

curves obtained from long-slit and two-dimensional Ha ob- 
■ 1 1 ■ i ■ ■ 

servations (see, e.g.,jd c Blok and Bosnia 2002; Swatcrs et alJ 
l2003l:ISimon et all2003Ll2005h . We shall restrict our analysis 
to these newer datasets in what follows. 

We also note that, strictly speaking, the observational 
disagreement is with the fitting formulae used to parame- 
terize the structure of simulated CDM halos (usually the 
profile proposed by NFW), rather than with the structure 
of simulated halos themselves. Although the fitting formu- 
lae provide a simple and reasonably accurate description 
of the mass profile of CDM halos, the radial range over 
which they have been validated often does not coincide 
with the scales where the disagreement has been identi- 
fied. This is important, since small but significant deviations 
between the NFW profile and simulated halos have been 
reported as the mas s and spatial resolution of the simula- 
tions has increased jMoore et al]ll99St iGhTgna et al-lliooot 
iFukushiee and Makinolll997l l200lF ~ In particular, the cen- 
tral slope of the density profile has been the subject of much 
debate. 

The NFW density profile is given by p(r) oc 
(f/f s ) _1 (l + r/r s )~ 2 , where r s is a scale radius. This pro- 
file ap proaches a power-l aw p (r) oc r" 1 "cusp" n ear the 
center. lMoore et all dl998f) and lGhigna et alJ ll200Cf) argued 
that halos simulated with higher resolution ha ve steeper cen- 
tral d ensit y profiles, p(r) oc r~ 1,5 . However. lHavashi et alJ 
l|2004h and lNavarro et alJ i2004h show that the density pro- 
file of halos becomes shallower towards the center and in 
fact may not actually converge to a well-defined asymptotic 
inner power law. At present, the simulation results of all 
groups are in good agreement, and the logarithmic slope 
of the density profile is ~ —1.2 or shallow er at the inner- 
most resolved radius JPiemand et alJl2004fl . Although there 
is no broad consensus yet regarding how deviations in sim- 
ulated halo profiles from fitting formulae like NFW may af- 
fect the comparison with observed rotation curves (see, e.g., 
lHavashi et al.ll2004l) . the fact that the deviations worsen to- 
wards the centre, in addition to the observed scatter in halo 
profiles, advise against using extrapolations of simple fitting 
formulae such as the NFW profile to assess consistency with 
observation. 

Finally, it must be emphasized that the "cusp vs. core" 
problem arises when comparing rotation speeds of gas in 
LSB disks to spherically-averaged circular velocities of dark 
matter halos. Given that CDM halos are expected to be sig- 
nificantly non-spherical, some differences between the two 
are to be expected. Previous studies have shown that sim- 
ulated halos are t riaxial objects wi t h shapes ranging from 
oblate to prolate jDj|vis_j3t_al||l985|; |Barnj3g_jaigLE[^ 
1987tlFrenk et al.ll988l:IWarren et al.ll992l : ljing et alll995t 
Thomas et al.lll998t IJing and Sutoll2002ri . The angular mo- 
mentum of the halo tends to be aligned with the minor axis, 
therefore one expects galactic disks to lie in a plane close 
to the plane defined by the intermediate and major axes of 
the halo. It is therefore important to take into account de- 
viations from spherical symmetry in the structure of CDM 



halos in order to make predictions regarding the rotation 
curves of gaseous disks that may be compared directly to 
observation. In fact, asymmet ries in the velocity fi elds of 
gas disks have been observed llCourteau et alj|2003l e.g.,), 
although in normal galaxies they are usually attributed to 
non-circular motion induced by bars or spiral arms. 

We address this issue here by exploring the closed loop 
orbits within simple non-axisymmetric potentials whose 
mass profiles are consistent with cold dark matter halos. 
We fo llow the formalism developed by Gerhard and Vietril 
( 1986) and iBinnev and Tremaind (11987 ) in the context o f 
tria xial bulges and bar potentia ls, and bv lFranx et all <ll994f) 
and ISchoenmakers et al.l <ll997l) to probe non-axisymmetry 
in spiral galaxies. We focus our analysis on the shape and 
inner slope of the rotation curves inferred for disks in triax- 
ial halos by simulating long-slit observations of their veloc- 
ity fields. We find that velocities along the slit may deviate 
significantly from the circular velocity, a result that may ac- 
count for the variety of LSB rotation curve shapes and that 
may reconcile the "cores" inferred for some LSBs with the 
cuspy mass profile of CDM halos. 



2 LSB ROTATION CURVES 

Figure Q illustrates the "cusp vs core" disagreement alluded 
to above. This figure shows the Ha rotation curves of two 
LSB g alaxies (points with e r ror ba rs) selected from the sam- 
ple of Ide Blok and Bosmal <200l B02) and ISwaters et all 
( 2003, S03). Also shown are the spherically-averaged V c pro- 
files of all galaxy-siz ed co ld dark matter halos presented in 
lHavashi et all &2004) and lNavarro et alJ (|2004T> . In order to 
emphasize discrepancies in shape, the rotation curves in Fig- 
urc^have been scaled to the radius, ro.3, and velocity, Vb.3, 
where the logarithmic slope of the curve is dlog V/d log r = 
0.3. The choice of this scaling is motivated by the asymp- 
totic slope of the NFW V c profile, d log V/d log r — 0.5. In 
this case, the scaling radius ro.3 = 0.411 r s is in the region 
well resolved by numerical simulation; it is also easily iden- 
tified in observed rotation curves, as it lies halfway between 
the linear (dlog V/d log r = 1) "core" region and the outer, 
flat (dlog V/d log r — 0) rotation curve portion seen in most 
disk galaxy rotation curves. 

The dashed line (which runs through the simulation 
data) shows the V c profile of an NFW halo, which is fixed 
and independent of concentration in these scaled units. We 
find that the shapes of the dark halo V c curves are in fact 
quite similar to NFW, although they span a large range of 
halo masses and have been chosen without regard to the 
dynamical state of the halo. 

LSB rotation curves, on the other hand, exhibit a wide 
range of shapes, and the two galaxies shown in Figure 
have been chosen to illustrate some of the extreme cases 
in the B02 and S03 datasets. These have been scaled to 
ro.3 and V0.3 as determined by fit s wit h the rotation curve 
fitting formula of lCourteaul (Il997f) (see lHavashi et aljf2004l . 
for details). Galaxies like UGC 5721 are roughly consistent 
with NFW and feature a gentle turnover between the rising 
and flat parts of their velocity curves. Others are better 
approximated by a linear scaling of velocity with radius, as 
expected in the presence of a constant density core. One 
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Figur e 1. Rotation curves o f two LSB galaxies from the sa mples 
o flde Blok and Bosmal <200l B02) and lSwaters et alJfcOOl S03), 
chosen to illustrate their various shapes. Rotation curves have 
been scaled to the radius 7*0.3 an d corresponding velocity Vb.3 
(see text for details). The NFW profile (dashed line) and the V c 
profiles of simulated dwarf- (dotted lines) and galaxy-sized halos 
(solid lines) match reasonably well systems like UGC 572f but 
cannot account for those with a linear rise in velocity with radius, 
like UGC 5750. 

such example, UGC 5750, is shown in Figure^ and is well- 
matched by the dot-dashed V oc r line. 



3 CLOSED LOOP ORBITS 

Does this result imply that a cusp in the dark matter density 
profile cannot be present in galaxies such as UGC 5750? As 
noted in §Q before concluding so one must take into account 
possible systematic differences between rotation speed and 
circular velocity in gaseous disks embedded within realistic, 
triaxial halos. This is a complex issue that involves a number 
of parameters, such as the degree of triaxiality, the role of 
the disk's self-gravity, size, and orientation, as well as the 
possibility of transient deviations from equilibrium. 

We explore a simplified scenario here, where the orien- 
tation of the principal axes of the potential is constant with 
radius and a massless, filled gaseous disk is placed on one of 
the symmetry planes of the halo. The velocity field of the 
disk may then be examined by solving for the closed loop or- 
bits within a two-dimensional non-axisymmetric potential. 
This is analogous to orbits within a barred potential where 
the pattern speed of the bar is set to zero, and has been 
studied by a number of authors using the epicycle approx- 
imation for the case of weak bars . Follo wing the formalism 
presented in ISchoenmakers et al.l lll997l) . we may write the 
potential as, 

$(7?, cj>) = <5> (R) + ® m (R) cos(roc/>o - <t>m(R)), (1) 
where &o(R) is the unperturbed potential, <& m (R) is 



a stationary perturbation to that potential, <f>o is the az- 
imuthal angle, and 4> m {R) is the radially-dependent phase 
of the perturbation. For the simple case we consider here, 
we set m = 2 and require the phase of the (stationary) 
perturbation to be constant with radius, cj> m {R) = <f> m . We 
then choose the origin so that 4> m = along the major (mi- 
nor) axis of the perturbed potential for positive (negative) 

The closed-loop orbit solutions can be calculated ana- 
lytically using the epicyclic approximation and are given in 
parametric form by 

-Rorb = Ro(l- ■■^p cos(m0o)) (2) 

, , dim T" &3m . , , \ (c .\ 

<p rb = <po H ~ sm(m0 o ) (3) 

2m 

V R = mV c {R ) ^ sm{m4> ) (4) 

U = V^ORo) (l + ^ cos(m0 o )) , (5) 

where (Ro,4>o) is the guiding centre of the closed loop 
orbit; <j>o = fiot; and the (unperturbed) angular and circular 
velocities are given by n% = 1/Ro d$/dR and by V C (R ) = 
«o.Ro. 

3.1 Case I: A Uniform Perturbation 

In general, the terms a\ m and az m 1 are functions of 4>o(-R), 
& m (R) and their first and second derivatives with respect 
to R. However, we begin by examining the simpler case of 
a uniform perturbation, i.e., $ m (i?) = f<&o(R), where / is 
a constant. In this case, these terms in the closed-loop orbit 
solutions simplify to: 
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(6) 



g(R) = $o7"o-l, (8) 

where the escape velocity is V CS c = y2|$o|; = 
d 2 &o/dR 2 ; and all terms in eqs. l6l8l are evaluated at R — Rq. 
We note that this perturbation tends to a finite value as 
R — > and therefore can result in unphysical solutions to 
Poisson's equation for the density distribution at small radii. 
However, we explore this simplified case in order to illustrate 
the general character of the closed orbit solutions. A more 
realistic perturbation in which / tends to zero as R — *• is 
presented later in this section. 

For a power-law potential, V c oc R a , the function g is 
independent of radius, g = 2(q — 1). Near the center of an 
NFW profile a ~ 1/2, and the ellipticity of the orbits scales 
simply as e R = (7? max - Rmm)/(Rmax + Rmitx) = Oim/2 = 
/ (Ksc/K 2 — 1). Note that the shape of the orbit at some 
radius depends on the mass profile both inside (through V c ) 
and outside (through 14 sc ) the orbit, since the gravitational 
forces of the outlying material do not cancel out in a non- 
spherical potential. Quantitatively, the deviation from circu- 
larity depends on the competition between the monopole of 

1 12m and ci4 m are zero for stationary perturbations 
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the unperturbed potential and the quadrupole of the pertur- 
bation. The larger the central mass concentration the larger 
the monopole and the smaller the ellipticity of the orbit 
for given fractional perturbation /. Indeed, deviations from 
circular motion are minimum for a point-mass (Kepler) po- 
tential, where er = /. For more general mass profiles, such 
as those of CDM halos, where V csc is finite at the center and 
V c tends to zero, the ellipticity of closed loop orbits is ex- 
pected to increase dramatically toward the center, even for 
relatively minor deviations from spherical symmetry. 

Along an orbit, the maximum deviations from circular 
motion occurs along the minor and major axis of the per- 
turbation. Assuming an NFW profile for the unperturbed 
potential, i.e., 3>o = $nfw, we consider in Figure the ra- 
dial dependence of the tangential velocity, V4,, for various 
values of /. 

The top panel shows the logarithmic slope of the V4, ve- 
locity profile along the (j>o — axis for m = 2 perturbation 
magnitudes of 0.05% and 1%, respectively. For negative val- 
ues of /, the (f)a = axis corresponds to the minor axis of the 
isopotential contours and the major axis of the closed loop 
orbits, and the tangential velocity is at a minimum at this 
point on the orbit. The middle panel of Figure [5] shows the 
ratio V<j,/V c , and the bottom panel shows the ellipticity of 
the orbits, which scales as the ratio between the escape and 
circular velocities and therefore increases toward the center 
of the NFW potential. 

Also shown in the top panel are the logarithmic slopes 
of the NFW (dashed) and pseudo-isothermal (solid) veloc- 
ity profile s , as w ell as the rotation curve fitting formula of 
ICourtearJ Jl997f) (dotted). The pseudo- isothermal velocity 
profile is given by V^ a (r) = V£(l - (r c /r) tan -1 (r/r c )), 
where Voo is the asymptotic velocity and r c is the radius 
of the constant density core in this model. The Courteau 
velocity profile is given by V 7 (r) = Vb(l + (r/Vt) -7 ) -1 / 7 , 
where Vb and r t are dimensional scaling parameters and 7 
is a dimensionless parameter that characterizes the shape of 
the rotation curve. This three-parameter formula provides 
excellen t fits to LSB rotatio n curves with a wide variety of 
shapes llHavashi et al.lEoo4) . 

Figure shows that even small perturbations in the 
potential can result in large changes in the shape of the 
velocity profile. For / = —0.005, the shape of the V$ profile is 
similar to that of the pseudo-isothermal profile (and of that 
corresponding to 7 = 1.5) down to radii as small as one-third 
of ro.3. Reducing the magnitude of / has the effect of moving 
inwards the radius (in units of ro.3) where the deviations 
become important but, because the Vg SC /Vc term formally 
diverges at the center, the shape of the rotation curve always 
deviates significantly from the true circular velocity curve 
somewhere near the center. Also shown in Figure [5] is the 
velocity profile along the (j>o = axis for a perturbation 
with / = 0.0005. Note that this is equivalent to the velocity 
profile along the minor (<f>o = 90°) axis for / = —0.0005 (see 
eqs. 5 and 7). Along this axis, the tangential velocity is at 
a maximum, and the logarithmic slope of the inner velocity 
profile is actually shallower than NFW. This implies that 
the shape of the rotation curve of a disk in a non-spherical 
potential varies depending on the viewing angle from which 
it is observed, a possibility we investigate further in0] 

We note that the analytic expressions for the closed or- 
bit solutions are derived assuming d(j> OI b/dt ~ dcj>o/dt, there- 
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Figure 2. Top panel: Logarithmic slope, dlogV/dlogr, of the 
pseudo- isothermal (solid ) and NFW (dashed) circular velocity 
profiles. Courtem]| ll997t) velocity profiles for three values of the 
parameter 7 (increasing from left to right) are shown by the dot- 
ted curves. All curves are normalized to the radius ro.3 where 
dlogV/dlogr = 0.3. Long-dashed, / < 0, (dot-dashed, / > 0) 
curves correspond to the tangential velocity profile, Va, measured 
along the major (minor) axis of the closed-loop orbits of a disk 
in an NFW potential with an m = 2 perturbation. The magni- 
tude of the perturbation is a constant fraction, /, of the NFW 
potential at all radii. Curves for two perturbation magnitudes are 
shown, with |/| increasing from left to right for the long-dashed 
curves. Even for rather small perturbations to an NFW potential, 
the tangential velocity profile exhibits a wide variety of shapes, 
matching the pseudo-isothermal profile over a wide range of ra- 
dius, as well as the Courteau profile with various values of the 
shape parameter 7. Middle panel: Ratio of to the NFW circu- 
lar velocity V c for the same two values of / as in the top panel. The 
departure from circularity depends on the ratio between escape 
and circular velocities (shown by the upper curve) The escape ve- 
locity becomes much larger than the circular velocity towards the 
center of the NFW potential, resulting in large deviations from 
circularity even for small perturbations in the potential. Bottom 
panel: Ellipticity of the orbits, e^j, for constant perturbation am- 
plitudes shown in the top panel. The ellipticity is simply related 
to the ratio between escape and circular velocities and therefore 
increases toward the center of the NFW potential. Note that the 
epicyclic approximation breaks down at small radii for sufficiently 
large perturbation magnitudes, resulting in invalid solutions with 
e R > 1. 



fore the solutions break down for large perturbations where 
epicyclic theory is not a good approximation to the true so- 
lution for the orbits. The bottom panel of Figure [2] shows 
that the ellipticity of the orbits predicted by the epicyclic 
approximation reaches values of ejj > 1 at small radii for 
sufficiently large perturbation magnitudes. In practice the 
closed orbit solutions are valid down to some minimum ra- 
dius where the epicyclic approximation still holds, therefore 
we show the solutions only at radii where d(j> or h / 'dt<,d(po / dt 
in all subsequent figures. 
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The results shown in Fig [5] suggest that caution should 
be exercised when interpreting the constraints on the slope 
of the dark halo inner cusp derived from analyzing the inner- 
most region of long-s lit rotation curves dde Blok and Bosmal 
120021: Ide Blokl l2005). These regions a re not just particu- 
larly sensitive to obs ervational error l|Swaters et alJ 2003; 
ISpekkens etafl l2005Tl but, as the discussion above illus- 
trates, one would expect on rather general grounds that gas 
velocities near the center might deviate substantially from 
the circular speed. 



3.2 Case II: A Radially Varying Perturbation 

In general, the magnitude of the perturbation, /, is expected 
to be a function of radius and it is instructive to consider 
perturbations to an NFW potential that may result in a 
pseudo-isothermal rotation curve along the major axis of 
the closed loop orbits. In order to derive such a perturba- 
tion, wc cannot use the simplified solutions in eqs. 16181 but 
instead set V0 = VS so and calculate numerically the function 
fiso(R) = $m/$NFW required to satisfy eq. ^Jfor cj>o = 0. 
Note that, by construction, this perturbation vanishes as 
i? — > since both the pseudo-isothermal and NFW mass 
profiles are convergent at small radii. 

We set the parameters of the unperturbed NFW po- 
tential to match a typical LSB galaxy; the mass, virial ra- 
dius and scale radius are set to M200 = 10 11 Ii~ 1 Mq, 
^200 = 75 /i _1 kpc, and r s = 6.25 /i _1 kpc, respectively, giving 
a peak circular velocity of V m ax = 95.2 km s" 1 . The pseudo- 
isothermal velocity profile is normalized to match the NFW 
profile at (ro.3> V0.3) = (2.57 /i -1 kpc, 73.5 km s _1 ) such that 
r c = 0.95 /i _1 kpc and Vac = 99.2 km s _1 . The maximum 
deviation between the NFW and the pseudo-isothermal pro- 
file within r < 20 /i _1 kpc is 11.7 km s _1 at a radius of 
r = 0.22 h~ kpc. We emphasize that, although this proce- 
dure specifies a particular radial dependence for f(R), in 
general most well-behaved perturbations will lead to de- 
viations from the circular speed which become more pro- 
nounced toward the center, generally mimicking the pres- 
ence of a "core" along the <j>o = axis. 

A few closed loop orbits are shown in the top panel of 
Figure [3] along with the isopotential contours of the per- 
turbed NFW potential. Also shown are the corresponding 
isodensity contours, calculated by solving Poisson's equa- 
tion. The major axis of the halo isopotential contours is per- 
pendicular to the major axis of the closed loop orbits. The 
bottom panel shows the eccentricity, e = 1 — Rmin/Rmaxt 
of the orbits as well as that of the isopotential and isoden- 
sity contours. In general, the isopotential contours are more 
circular than either the isodensity contours or the orbits, 
and all three become more elongated at small radii, despite 
the fact that the magnitude of the perturbation decreases 
towards the center. 

This is shown in the bottom panel of Figure which 
shows the perturbation magnitude as function of radius. The 
magnitude of the perturbation fi so (R), increases from zero 
at the center to a maximum of 0.375% at r ~ 0.25 ro.3 and 
is well fit by a function of the form: 

f(R) = axexp(-x/b), (9) 
where x = R/r 3 and the best fit parameter values are 
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Figure 3. Top panel: Closed loop orbit trajectories in a perturbed 
NFW potential, as well as isopotential and isodensity contours. 
Orbits are elongated in the direction perpendicular to the major 
axis of the isopotential contours. Middle panel: Eccentricity of 
the orbits, as well as that of isopotential and isodensity contours, 
as a function of the guiding centre radius, Ro. Bottom panel: 
Perturbation profile, /i so (solid curve), and fitting formula given 
by eq. [5] (dashed curve) . 



a = 0.1 and b = 0.098. We note that the perturbation de- 
scribed by eq. |3] or any similar function which peaks at 
roughly the same radius and amplitude, gives similar results 
as the exact solution for fi so (R). 



4 SIMULATED ROTATION CURVES 

We now investigate the rotation curves expected when ob- 
serving a thin gaseous disk whose kinematics are described 
by the closed loop orbits presented in the previous section. 
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Figure 0| shows the orbits for a disk projected at an inclina- 
tion angle i = 60° . The upper panels of each of the four sets 
show the projected image of the disk, color-coded by line- 
of-sight velocity Vi os . The isovelocity contours are drawn at 
10 km s~ intervals. The lower panels show the rotation 
curve obtained by sampling the disk with a thin slit placed 
along the major axis of the disk, as indicated in the up- 
per panels by the solid horizontal line. Also shown are the 
NFW and pseudo-isothermal velocity profiles and rotation 
curve data from one of the two galaxies in Figure scaled 
to (ro.3, Vo.3). 

The upper left, upper right, and lower left sets of panels 
of Figure0]show, for / = fi so (R), the disk when viewed along 
various orientations, so that the angle between the projected 
minor axis and the minor axis of the orbits is cj) a = 0° , 15° , 
and 30°, respectively. 

When the two minor axes coincide, V\ os (R) = 
Vtp(R,(f)o — 0°)sin(i), and the rotation curve matches the 
Viso curve by construction (see open circles in the rotation 
curve panel). Note that for this projection, the isovelocity 
contours of the disk are exactly symmetric, and the ellip- 
tical nature of the orbits is degenerate with the inclination 
and circular velocity profile, so that the triaxiality of the 
halo would be difficult to detect. Indeed, a slit placed across 
the minor axis of the disk — a crude but useful method of 
assessing triaxiality — would show no net rotation for this 
projection (see diamonds in the rotation curve panel). This 
ought to be taken into account when searching for signa- 
tures of tria xiality in 2D velocity fields, and it weakens the 
argument of lGentile et alJ 1120051) . who argue against triax- 
iality from the lack of strong minor axis rotation in the HI 
velocity field of DDO 47. Actually, the isovelocity contours 
in DDO 47 (see their Figure 1) show clearly some of the 
twists and asymmetries expected for elliptical orbits, as we 
discuss below. 

Under a different projection, (f> a — 15°, subtle differ- 
ences in the velocity field are noticeable. Firstly, the isove- 
locity contours are no longer symmetric about the x axis. In- 
stead they exhibit "twists" near the centre of the projected 
disks, although the velocity field is symmetric relative to a 
180° rotation. As a result, the minor axis no longer traces 
an isovelocity contour, and some "rotation" is seen along 
the minor axis, although confined to the very center. Note 
that the shape of the major axis rotation curve still matches 
quite well the pseudo-isothermal profile and the linearly ris- 
ing rotation curve of UGC 5750. 

The lower left panels of Figure 0] shows the disk from 
a more skewed viewing angle, <j} a = 30°. The major axis 
rotation curve in this case resembles better the true NFW 
circular velocity profile and the rotation curve of UGC 5721. 
The distortions in the isovelocity contours are now more pro- 
nounced; in particular, the shape of each contour is different 
in the upper half from the bottom half of the galaxy but it is 
still symmetric relative to a 180° rotation. Loosely speaking, 
one may see the imprint of the halo triaxiality in a velocity 
field that differs significantly in contiguous quadrants, but 
is anti-symmetric in diagonally opposite ones. 

Complex velocity fields requiring more than just circu- 
lar, coplanar motion in "tilted rings" are the norm rather 
than the exception near the center of LSBs. This has 
been clearly demonstrated by the high-resolution combined 
HI, Ha, and CO maps of 5 LSBs presented recently by 



ISimon et alJ J2003L 120051. As these authors conclude, depar- 
tures from simple circular motion is clear in the majority of 
galaxies in their sample. 

Although Simon et al interpret these deviations in term 
of "radial motions" (i.e., pure expansion or contraction com- 
ponent of the tilted rings) the interpretation of this complex- 
ity in terms, instead, of elongations in the potential remains 
viable. In particular, we note that the isovelocity contours 
of NGC 2976 (Figure 4a of Simon et al 2003) show clearly 
some of the twists illustrated in the lower left panel of Fig- 
urc0| We also note that their conclusion about the lack of a 
cuspy dark matter halo relies heavily on assuming that the 
rotational velocity component of their tilted-ring fit traces 
the circular velocity of the halo, an assumption that may re- 
quire revision if the non-circular motions seen in their data 
are ascribed to the effects of triaxiality. Ongoing work is 
aimed at reaching a more definitive conclusion on this. 

Finally, the lower right panels of Figure[I]show, for com- 
parison, the velocity field of a disk in circular motion in an 
unperturbed (/ = 0) NFW potential, whose rotation curve 
is in this case indistinguishable from the long-slit, / = /i so , 
(f> — 30° rotation curve. This illustrates that long-slit ro- 
tation curves cannot discriminate between these two cases; 
two-dimensional velocity fields are needed. 



5 DISCUSSION 

Given the variety of rotation curve shapes that arise when 
considering gaseous disks in non-axisymmetric halo poten- 
tials, the outlook for reconciling dark matter cusps with LSB 
rotation curves is promising. In particular, the fact that de- 
partures from circularity are expected to increase towards 
the center, and that relatively large effects are possible even 
in the presence of minor asphericity appear rather encour- 
aging. 

However, it would be premature to argue that the prob- 
lem has been fully solved. Although from the analysis pre- 
sented here it appears that, on general grounds, one may 
expect at least some rotation curves to be consistent with 
cored profiles, it is unclear at this point whether the fre- 
quency (and sizes) of apparent "cores" is consistent with 
this interpretation. This is likely to depend critically on the 
triaxiality of the potential in the very inner regions, as well 
as on the response of the halo to the assembly of the disk 
and its orientation relative to the halo's principal axes. All 
of these issues are still quite uncertain. 

It is therefore important to build a more compelling 
case for the triaxial interpretation of cores in LSB rotation 
curves, so as to render it falsifiable. Are the shapes of CDM 
halos in fact consistent with the perturbed potentials we 
have explored? We are currently investigating this using a 
sample of high-resolution galaxy-sized halos. Are there any 
corroborating traits that may be used to confirm or exclude 
the hypothesis that halos surrounding LSBs are indeed tri- 
axial? Baryonic matter in the form of a bar could have a 
simil ar effect on g as motions as a triaxial dark matter halo 
jRhee et al.ll2004l1 . 

How can one best verify the triaxial halo interpreta- 
tion in two-dimensional velocity maps? Identifying a clean 
and unambiguous indication of triaxiality, such as the un- 
usual isovelocity "twists" and symmetry patterns shown in 
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Figure [I] will be as important as the success of aspherical 
halos in reproducing the rich variety of shapes of LSB ro- 
tation curves. Only if this is accomplished shall we be able 
to conclude that LSB rotation curves do not preclude the 
presence of dark matter density cusps, thereby freeing the 
CDM paradigm of one vexing challenge on small scales. 
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Figure 4. Upper panels of each set: Projected velocity fields of disks in perturbed (/ = /i so ) and unperturbed (/ = 0) NFW potentials, 
inclined by i = 60°. The disks are color-coded by line-of-sight velocity as indicated by the color bar. The green curves show isovelocity 
contours in 10 km/s intervals. The solid black line indicates the position of a long slit placed along the major axis of the projected disk. 
The (/ = /i so ) disk is shown for three different values of <j> a , the angle between the slit and the long axis of the closed loop orbits. Lower 
panels of each set: Rotation curves corresponding to the disks in the upper panels, scaled to (rrj.3, Vo.3)- Open circles (open diamonds) 
show the velocity profile along the slit placed across the major (minor) photometric axis of the projected disk. Open triangles (squares) 
show the rotation curve of UGC 5750 (UGC 5721) from Figure^ The shape of the rotation curve of the / = /i so disk agrees with that 
of UGC 5750, which rises linearly with radius for <f> a = 0° and 15°, and is similar to the pseudo-isothermal velocity curve shown by the 
solid curve. The minor axis velocities are zero when the disk is viewed at (fr a = 0, but nonzero minor axis velocities and twisted isovelocity 
contours result when a different line-of-sight, (fr a = 15° or 30° is chosen. The / = /i so , <t> = 30° and / = disks have major axis rotation 
curve shapes similar to or identical to the NFW velocity profile, and match closely the UGC 5721 rotation curve. The velocity field of the 
disk in the perturbed potential can be distinguished from the unperturbed NFW by the nonzero minor axis velocities and the twisted 
isovelocity contours. 



